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ABSTRACT 

The  general  mathematical  problem  of  the  thermal  entrance 
region  Is  formulated  for  a  parallel  plate  channel  by  Including  the 
effects  of  viscous  dissipation,  axial  conduction,  and  wall  blowing. 

The  associated  eigenvalue  problem  Is  solved  by  the  B.  G.  Galerkin 
method  and  the  results  are  presented  for  constant  wall  temperature. 

It  is  shown  that  the  particular  method  has  distinct  computational 
advantages  over  the  classical  form  of  solutions.  The  constant  wall 
temperature  case  Is  investigated  by  employing  the  solutions  of  the 
eigenvalue  problem,  and  It  Is  concluded  that  the  wall  blowing  param¬ 
eter  has  considerable  effect  on  the  temperature  development.  The 
axial  conduction  term  Is  also  shown  to  have  considerable  effect  on  the 
temperature  development  for  low  values  of  Peclet  number. 
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CHAPTER  I 
INTRODUCTION 


The  heat  transfer  In  laminar  flow  of  an  Incompressible,  non- 
viscous  fluid  in  the  entrance  region  of  a  circular  tube  was  first 
Investigated  analytically  by  Graetz  (I)1,  under  the  assumption  of 
specified  uniform  wall  temperature  and  fully  developed  velocity  pro¬ 
files.  Prlns,  Mulder,  and  Schenk  (2)  solved  the  analogous  problem  with 
plane  parallel  plates  substituted  for  the  cylindrical  tube  wall. 

Various  authors  (3-8)  have  extended  analysis  of  the  problem  by  considering 
jincreaslng  complexity  of  the  boundary  conditions  and  velocity  distri¬ 
butions  for  circular  tubes  and  parallel  plates. 

Abramowitz  (9)  and  Sellars,  Tribus,  and  Klein  (10)  have  pre¬ 
sented  In  detail  the  classical  mathematical  solutions  of  the  eigenvalue 
problem  encountered  In  the  analysis  of  thermal  entrance  regions.  These 
authors  concluded  that  the  evaluation  of  the  eigenvalues  and  the 
eigenfunctions  from  the  exact  analytical  expressions  resulted  in  con¬ 
siderable  computational  difficulties  even  in  the  cases  of  relatively 
simple  flow  and  thermal  boundary  conditions. 

Sparrow  and  Siegel  (II)  introduced  the  successful  application 
of  variational  methods  to  the  mathematical  problem  for  parallel  plate 
channels  and  finite  cross-sectional  ducts  and  approximate  solutions 

^Numbers  in  parentheses  refer  to  similarly  numbered  references 
in  the  bib  I iography. 
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were  found  to  be  In  reasonable  agreement  with  the  available  exact  solu¬ 
tions  for  certain  special  cases.  Direct  numerical  solutions  for  various 
hydrodynamic  velocity  profiles  in  circular  tubes  were  presented  by 
Kays  (12)  who  compared  the  results  for  the  Graetz  problem  against  the 
exact  solutions. 

Hsu  (13)  presented  an  exact  mathematical  solution  for  the 
entrance  region  laminar  heat  transfer  Including  the  effects  of  axial 
conduction.  The  temperature  solution  for  pipe  flow  corresponding  to  the 
boundary  condition  of  uniform  wail  heat  flux  was  obtained.  The  corre¬ 
sponding  problem  in  parallel  plate  geometry  was  treated  by  Agrawal  (14) 
for  the  boundary  conditions  of  constant  wall  temperature.  The  eigen¬ 
functions  were  expanded  as  an  infinite  Fourler-slne  series,  and  the 
temperature  solution  corresponding  to  a  Peclet  number,  Pe  =  I ,  based 
on  the  mean  velocity  was  obtained.  These  authors  concluded  that  the 
effect  of  axial  conduction  was  significant  for  low  Peclet  numbers. 

LeCroy  (15)  investigated  the  feasibility  of  application  of 
the  B.  G.  Galerkin  method  (16)  to  formulate  approximate  solutions  for 
the  temperature  development  in  parallel  plate  MHD  channel  flows  by  in¬ 
cluding  the  effects  of  viscous  dissipation.  Joule  heating,  non-zero  net 
current,  and  axial  heat  conduction. 

The  previously  mentioned  authors  (1-15)  Investigated  the 
thermal  entrance  region  problem  assuming  channels  and  ducts  with  solid 
walls. 

The  purpose  of  the  present  study  is  to  develop  approximate 
solutions  for  the  temperature  development  and  local  Nusselt  number  in 
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parallel  plate  channel  flows  with  wall  blowing  by  the  application  of  the 
B.  G.  Galerkln  method  (16).  The  problem  was  formulated  to  include  the 
effects  of  viscous  dissipation  and  axial  heat  conduction. 
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CHAPTER  II 

ANALYTIC  PROCEDURE 

I.  GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 

Consider  a  fluid  of  constant  properties  flowing  between 
infinite  parallel  plates.  The  mathematical  problem  Is  formulated  In  the 
Cartesian  coordinate  system  (X|,  Xj,  x^),  with  x3  taken  at  the  base  of 
the  channel  along  the  direction  of  the  applied  pressure  gradient, 
dP/dxj.  The  origin  of  the  system  Is  taken  at  the  point  of  Intersection 
of  the  Xj  axis  with  the  plane  at  which  the  thermal  entrance  region 
Initiates,  as  shown  In  Figure  I.  It  is  assumed  that  there  is  a  constant 
injection  of  fluid  through  the  wall  at  X2  =  0,  and  that  there  Is  a  con¬ 
stant  suction  of  fluid  through  the  wall  at  X2  =  L.  The  mass  flow  rate 
due  to  Injection  Is  assumed  to  be  equal  to  that  due  to  suction  producing 
a  velocity  component,  Vw.  It  is  further  assumed  that  the  flow  Is  fully 
developed  at  the  thermal  entrance  plane.  • 

For  fully  developed  flow  conditions  and  for  a  parallel  plate 
channel  configuration  of  Infinite  dimensions  along  X|  and  x3,  with  no 
velocity  component  along  Xj  direction.  It  Is  concluded  that  the  applied 
axial  pressure  gradient,  dP/dx3,  is  a  constant. 

The  governing  equations  for  fully  developed,  steady  state. 
Incompressible  flow  of  a  viscous  fluid  with  constant  properties  for 
Infinite  parallel  plates  along  the  Xj  direction  are: 
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Continuity  equation, 


KJ  <x2>  *  0 


which  integrates  into. 


V2<x2>  =  Vw  =  constant 


Momentum  equation, 


32V3 

3X? 


(x2> 


3V3 

Vw  ^  (X2) 


3P 


pSx-, 


Energy  equation. 


P  cP 


Vw  3j<2  (x2*x3)  +  V3(x2)  3x^  (x2*x3> 


32T 


3^T 


3,x.  (x2»x3>  +  ^7  (x2>x3> 

2  3 


+  9  (x2> 


where  the  dissipation  function,  t(x2),  is  given  by. 


(I) 


(2) 


(3) 


(4) 


*(X2) 


(5) 


The  hydrodynamic  boundary  conditions  on  the  velocity  component, 
^3^x2^*  are  determined  according  to  the  nonslip  conditions  on  the 
stationary  wails. 


V3<0)  =  0,  V3(L>  =  0 


(6) 
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The  temperature  distribution  Is  prescribed  at  the  thermal 


entrance  plane,  x3  =  0,  as  a  function  of  x2, 

:  si  stews'!  60,  bns  aeldeHev  aaelrtolansniib  pr.lwcllot  oiiT  itmw 

T(x2,0)  =  T0(x2)  (7) 

<S*\T_V  _  S*  _ 

and  Is  prescribed  as  equal  ~tj^pera+\jfe&  at5  bath”  wal  Is’  of -the  ^aral  lei 


plate  channel. 


wT  "  oT 


-  <£X<£*)T 


1  — ■y--"— ^TjC  0.  , 

T(0,x3)  =  'tw,~m,x3>  =  Tw 


Tf 
O ' 


-  n 


(8) 


The |  Remaining  condition  on  the  temperature  -^pw^^onfyls* determined 

from  the  thermally  fully  developed  condition,  T,,^, x3),  by  solving 

-•isgmst  eoneitne  silt  tli  iw  beteioosaa  eulsv  oitanetoenedo  s  ai  *7  a-ierlv. 

the  limiting  case,  x^  -*■  »,  of  the  energy  equation,  ° 

. (ox)0T  ^noltudlntatb  a-tuta 

,£  no  I  fsup3  ^ngii&p^miytiiomgnyent.^^  or!]'  jgj 

2  p  w  dx_  2  2 

emooid  * y  I av  I i'osc53»"i  r.oltsupJJ  ..no-taupe  yjnorio  silt  bne 

under  the  boundary  conditions. 


(£1) 


VO)  =  Tw,  TJL)  =  Tw 


,  ,  l>b  „  .  .  L>D 

=,CV>  yb  V’  "  ^  3^ 


(10) 


There  would  be  little  change.fn  tire  further  mathematical  development  of 

(s,y)  p.  (j.ylS*  *  <v;»  eS 

the  problem  if  the  thermal  bouridariy  conditions  on  the^channel1  walls  were 

assumed  to  be  constant  wall  temperatures  qt. unequal  value. 

<M  )  Cy >V  +  's4y)  + 

II.  NON-DIMENSIONAL  PROBLEM 

?l  noltonut  noltcqlceib  aaejnoianemtb  eriT  e-jsflw 

The  mathematical  system  Is  non-dimensional i zed  based  on  the 
cpigsjtjant  pressure  gradient  by  defining  ajq^r^lliefrllsIiliCBve^oeity,  Vc, 
as 
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V  s  _  L2  dP 
Vc  p v  ^ 

with  the  following  dtjnenslonless  variables  and  parameters: 


(II) 


*2  *3  V3(X2> 

y  ■  r~  #  z  ■  r-  •  u(yJ  =  — v —  * 


o  ..  T(X2,X3>  -  Tw  a  _  T0  ~  Tw 

O  —  —  ViL  «i  $  0, 


T*  -  T 

o  w 


0  “  T*  T  ' 

o  * 


*(yJ  =Rrr*~rw)  *(*2> 


(12) 


where  T*  Is  a  characteristic  value  associated  with  the  entrance  temper¬ 
ature  distribution,  T0(x2). 

The  dimensionless  forms  of  the  momentum  equation.  Equation  3, 
and  the  energy  equation,  Equation  4,  respectively,  become 


d2u 

dy7 


du 


(y)  -  Rc\  3y  (y)  a  -  I 


(13) 


Pe  [u<y>  ||  <y,z>  +  X  ||  (y,z)j, 


a2e 

17  (y'z) 


g2g 

+  — r  (y,z)  +  v(y) 
3z 


(14) 


where  the  d  Intension  jess  dissipation  function  Is 


¥(y) 


Pr  Ee  [37  «»>] 


(15) 
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In  Equations  13  through  15  the  dimensionless  parameters  are  defined  as: 

V  ^ 

Ec  =  ^  =  Pseudo.  Eckert  Number 

=  Pec  let  Number 

=  Prandtl  Number 

=  Wall  Velocity  Parameter 

=  Characteristic  Reynolds  Number 

=  Wall  Blowing  Parameter  (16) 

mens  ion  I  ess  velocity  u(y),  becomes 

u(0)  =  0,  u(l)  =  0  (17) 

and  the  thermal  boundary  conditions  for  the  dimensionless  temperature, 
0(y,z),  at  the  thermal  entrance,  and  at  the  walls  are: 

0(y,O)  =  0o(y),  0(0, z)  =  0,  0(L,z)  =  0  (18) 

The  remaining  boundary  condition  required  to  complete  the  solu¬ 
tion  is  determined  from  the  thermally  fully  developed  dimensionless 
temperature  distribution,  0OB(y), 

0„(y)  =  IIm  0(y,z)  (19) 

2  -*■  m 


r  o 


Pe  = 


PVc^pL 


Pr  = 


Cp  pv 


A  =  JL 


Rc  = 


Vc  L 


RcX  = 


Vw  L 


The  boundary  condition  for  the  di 
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Solving  the  limiting  case,  z  +  •,  of  the  non-dimensfona!  energy  equation. 
Equation  14,  the  mathematical  system  becomes, 

d2e«  de^ 

- 5-  (y)  -  Pr  RcA  -7—  (y )  t  y(y)  =  0 

dy 

e»<0)  =  °»  eJ'>  =  0  (20) 

The  solution  for  the  fully  developed  dimensionless  temperature  distri¬ 
bution,  Equation  20,  Is  given  In  Appendix  B. 

III.  GENERAL  METHOD  OF  SOLUTION 


Since  the  fluid  properties  are  assumed  to  be  constant,  the 
dimensionless  momentum  equation  for  the  velocity  profile,  u(y),  is 
Independent  of  the  energy  equation.  The  solution  of  Equation  13  for 
u(y)  is  given  In  Appendix  A  as 


u(y) 


I  -  e^Y 
I  -  eRc* 


(21) 


Hence,  the  dimensionless  energy  equation  represents  a  linear,  non- 
homogeneous,  partial  differential  equation  with  variable  coefficients 
subject  to  the  given  boundary  conditions. 

Defining  a  new  independent  variable,  £,  as  the  dimensionless 
axial  coordinate  scaled  by  the  Pec  let  number  and  a  new  dependent 
variable,  fl(y,z),  as  the  dimensionless  excess  temperature  which  Is  the 
difference  between  the  local  temperature  and  fully  developed  temperature 
distribution,  6^,  as: 
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n(y,5)  -  6(y,z)  -  Ojy)  (22) 


the  energy  equation,  upon  substitution  of  Equation  22  Into  Equation  14, 
becomes 


u(y)  |£  (y,5)  +  PeX|^(y,5>  -  *2° 


i  a2ti 


H 

320. 

ay2" 


a  l  <y,5)  "  pT  — ?  (y»?) 
3y*  Pez  3^ 


30. 


ty)  -  PrRcx  ^y-  (y)  +  v(y) 


(23) 


The  bracketed  term  on  the  right  hand  side  of  Equation  23  Identically 
vanishes 'according  to  the  fully  developed  solution  of  the  mathematical 
system.  Equation  20,  hence  the  differential  equation  for  the  dimension¬ 
less  excess  temperature,  fi(y,£),  reduces  to  a  homogeneous  partial 
differential  equation  with  variable  coefficients. 


u(y)  (y,C)  t  PrRcX  <y,5)  =  (y,£)  +  >  (y,£)  (24) 

4  dy  3y^  Pe^  Hz 

with  boundary  conditions: 

ft(y,0)  =  ft0<y)  =  0o(y)  -  e^Cy) 

fl(0,E)  =  0,  0(l,C)  =  0 

l,m  =  o  (25) 

£  -*■  00 
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The  general  solution  of  the  differential  equation.  Equation 
24,  can  be  constructed  as  an  infinite  series  of  the  form: 


A  ,  .  -6nC 

fi(y,C)  =  V  a1"'  ftn(y)e 

n=*l 


(26) 


where  the  eigenfunctions,  ftn(y),  and  the  associated  eigenvalues,  0n, 
are  to  be  determined  from  the  homogeneous  eigenvalue  problem. 


(y)  -  PrRcA 


dnn  T  ^n2 

37-  (v>  +  [en  uty)  +  ^]Vy>  = 


with  the  boundary  conditions 


(27) 


nn(0)  =  o,  nnu)  =  o 


(28) 


A  simplification  of  Equation  26  is  possible  by  defining  a  new  set  of 
eigenfunctions  of  the  form, 


F  (y)  =  e“!sPrRcxY  £2  (y) 


n ' ' '  "  '  n 

Hence  the  eigenvalue  problem  reduces  fo. 


(29) 


d2F, 

dy' 


n 


(y)  + 


Bn2  Pr2Rc2A2 


Fn(Y)  =  0 


(30) 


Fp(0)  =  0,  Fp(l)  =  0 


The  general  solution  of  the  differential  equation.  Equation 
24,  becomes 
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JlCy,5) 


CO 


i 

n=  I 


a(n) 


Fn(y)e 


^PrRcXy 


e 


(31) 


Since  axial  conduction  and  wall  blowing  effects  are  Included, 
the  general  form  of  the  differential  equation.  Equation  30,  does  not 
satisfy  the  conditions  of  the  classical  Sturm-LIouvl I le  system  with  a 
complete  set  of  eigenfunctions  which  are  orthogonal  with  respect  to 
the  weighing  function  u(y).  If  we  assume  that  Fn(y)  and  Fm(y)  are  two 
characteristic  functions  which  are  solutions  of  Equation  30  and  let  Bn 
and  0m  be  the  corresponding  eigenvalues,  assuming  m  j*  n,  then  the 
integral  relation  between  the  associated  eigenfunctions  becomes 


I 

<s„  -  W  /  uty>  Fn<y>  Fra<y>  «y 

O 

Fn(y>  Fm(y)  dy  =  0  (32) 

o 

Equation  32  reduces  to  the  Strum-LIouvl I le  system  satisfying  the  pro¬ 
perty  of  orthogonality  for  the  limiting  case  of  Pe  ■+■  «,  hence  the  present 
problem  may  be  considered  as  a  more  general  class  of  eigenvalue  problem 
of  thermal  entrance  regions. 
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CHAPTER  III 

APPROXIMATE  SOLUTION  OF  THE  EIGENVALUE  PROBLEM 
BY  B.  G.  GALERKIN  METHOD 


The  analytical  construction  of  the  exact  solutions  for  the 
eigenvalue  problem  of  thermal  entrance  regions  is  difficult  even  In  the 
case  of  hydrodynamic  conditions  with  the  axial  conduction  term  and  wall 
blowing  being  neglected.  The  eigenfunctions,  for  this  restricted  class 
of  Sturm-Liouvi I le  systems,  are  obtained  as  infinite  series  of  Bessel's 
functions  of  various  orders  which  create  considerable  computational 
difficulties.  Since  the  present  problem  represents  a  more  generalized 
class  of  eigenvalue  system,  one  expects  at  least  the  same  amount  of 
mathematical  complexities  in  the  construction  of  the  exact  solution. 
Although  numerical  solutions  of  the  eigenfunction  differential  equation, 
Equation  30,  are  possible,  the  two  point  boundary  conditions,  together 
with  the  unknown  Bn,  necessitates  simultaneous  iterative  schemes  for  the 
eigenvalues  and  eigenfunctions  which  result  in  extensive  computational 
effort. 

It  is  proposed  to  obtain  a  finite  set  of  approximate  eigen¬ 
values  and  associated  eigenfunctions  for  Equation  30  by  the  application 
of  the  B.  G.  Galerkin  method.  Let  |'t|<(y)j  b®  a  selected  set  of  functions 
with  continuous  first  and  second  order  derivatives  in  the  range, 

0  S  y  S  that  satisfy  the  homogeneous  boundary  conditions  identically 
for  a:i  k.  The  selected  function  set  must  satisfy  the  conditions 


14 


AEDC-TR. 69-58 


♦k(0)  =0,  <fk(l)  =  0  (k  =  1,2 . N)  (33) 

The  approximate  solution  for  Fn (y )  is  constructed  as  a  linear 
combination  of  the  finite  set  of  functions  {$k(y)>  as: 

N 

Fn(y)  =  <j> j (y)  =  a(<n)  4> ( <y )  +  ...  ak(n)  ^(y)  (34) 

i  =  l 

Since  the  set  <|>k(y)  identically  satisfies  the  boundary  con¬ 
ditions,  Equation  33,  for  all  k,  C y )  also  satisfies  the  homogeneous 
boundary  conditions.  Substituting  Equation  34  Into  Equation  30  gives 
the  following. 


The  non-dimensional i zed  velocity  profile,  u(y),  is  not 
necessarily  symmetric  about  the  channel  centerline;  therefore,  the 
eigenfunctions,  Fn(y),  are  expected  to  yield  non-symmetrical  solutions. 
Considering  the  form  of  the  boundary  conditions.  Equation  33,  a 
trigonometric  function  set,  {^(y)},  is  selected  as: 

$k(y)  =  sin  k  tt 

k  =  1,2,3,  -  N  (36) 

Following  the  method  of  B.  G.  Galerkin,  that  Is,  multiplying 
Equation  35  by  the  functions  of  the  set,  <f>k£y)  =  sin  kiry,  and 
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integrating  over  the  range,  0  *  y  *  I,  the  linear  system  of  equations 
becomes: 


T. 

!  _  I 


'i(n)  Vi<M> 


Pr2Rc2x2 

— n — 


I2,2 


i2(k,n 


=  0 


(37) 


The  expressions  for  the  integrals  in  Equation  37  can  be  obtained  by 
analytical  Integration  using  the  velocity  distribution,  u(y),  given  in 
Equation  21  as: 


if  i  ¥  k. 


Ij(k,i)  =  J  u(y)sin  i Try  sin  kiry  dy  = 
o 


I  j  I-  <-l)i+k  f  I 


^  ( i+k)2  ~  ( i-k)2 


- - -  [l-  t-l)i+keRcAl  f  -- - ! 

2(eRcA-|)  L  J  [  Rc2X2  + 


(i+k)2ir2 


Rc2A2+(i-k)2,2 


if  i  =  k. 


l2<e« 


cX-l))['  " 


(eM  -|) 


RCA  R  2A2  +  4 1  2tt2 


(38) 


(39) 
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where  6jk  Is  the  Kronecker  delta.  Therefore,  Equation  37  becomes  a 
homogeneous  system  of  linear  equations  for  the  undetermined  linear 
combination  constants,  a}(n).  If  a  non-trivial  solution  is  to  exist 
for  the  aj(n*'s,  hence  giving  Fn(y)  a  non-trivial  solution,  the 
determinant  of  the  coefficient  matrix  must  vanish  identically  for  each 


Pr2Rc2X2 

4 


(41) 


Equation  4:  represents  a  polynomial  of  at  most  2N  degree  in 
0n,  with  possibly  N  distinct  positive  roots.  A  set  of  (N-l)  linearly 
independent  equations  for  a^n^'s  exists  for  each  evaluated  as  a  root 
of  Equation  41  which  may  be  solved  in  terms  of  one  of  the  unknown  linear 
combination  constants,  a^1'*,  whiqh  will  remain  undetermined.  Selecting 
the  particular  undetermined  constant  as  a^0^  =  An,  the  linear  system  of 
equations  for  the  remaining  constants  a.*n^  becomes: 
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where 


c,<">  .  fl 


(n) 


a,(n) 


al 


TnT  = 


k  =  2,3, ...N 


(43) 


Equation  42  is  a  non-homogeneous  system  of  linear  equations 
which  can  be  solved  to  determine  the  Cj^’s  for  each  $n  evaluated  as 
a  root  of  Equation  41.  Hence  the  approximate  eigenfunction,  Fn(y), 
associated  with  8  becomes 


Fn(y)  =  An  cj(n)  sin  Iiry 

i  =  l 


(44) 


where  cj^n^  =  I,  and  the  Gj^n*'s  are  known  from  the  solution  of  Equation 
42,  Substituting  Equation  44  Into  Equation  31  yields  the  approximate 
solution  for  the  dimensionless  excess  temperature,  fl(y,£),  as: 


K  f  N 

a<y,e>  =  £  Apje^Rc^  y  c,<")  sin 
n=l  L  1=1 


(45) 


where  K  is  the  maximum  number  of  distinct  eigenvalues  evaluated  as  the 
roots  of  the  determinant,  Equation  41. 

By  applying  the  boundary  conditions  given  in  Equation  25  the 
remaining  K  number  of  constants,  An,  can  be  determined  from 
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»(0,y) 


cj(n)  sin  Iiry 


=  n0(y) 


(46) 


which  represents  the  function  0Q(y)  as  a  linear  combination  of  known 
continuous  functions,  0n(y),  In  the  range  0  *  y  *  I.  Since  the  finite 
set  $j(y)  =  sin  iiry  satisfies  the  necessary  regularity  conditions,  An 
can  be  determined  by  the  application  of  an  extension  of  the  Weierstrass 
approximation  theorem  (17),  That  is,  multiplying  Equation  46  by  nm(y) 
and  integrating  over  0  *  y  ^  I,  the  system  of  equations  becomes: 


K 


T.  An 


CKm,n)U  =  L(m) 


(47) 


n=l 


where  the  integrals  are  defined  as: 


I  r  N  -|  r  Jt 

I(m,n)  =  J  ^  cj(n)  sin  iiry  2_, 


sin  jiry 


LJ-i 


PrR-Xy 

e  dy  (48) 


L(m) 


where 


1 

r 

r  n 

J  fl0(Y> 

I  =j(m) 

L  i  =  l 

sin  iny 

o 

-  c  (n) 
cl 

=  1;  a0(y> 

=  eo<Y> 

JsPrRcAy 
e  dy 


(49) 


(50) 


The  integrals  of  Equations  48  and  49  may  be  evaluated  by  analytical 
integration  using  the  fully  developed  temperature,  0OB(y),  as  presented 
in  Appendix  B, 
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Since  C|tm),s  were  previously  determined.  Equation  47  repre¬ 
sents  a  non-homogeneous  system  of  K  linear  equations  In  K  unknowns  which 
can  be  solved  for  the  An's  to  complete  the  approximate  solution  for  the 
excess  temperature  distribution. 
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The  general  method,  developed  In  the  previous  chapter,  was 
applied  to  the  associated  eigenvalue  problem,  Equation  30,  with  two 
finite  sets  of  trigonemetric  functions,  N  =  8  and  N  =  16.  Using  the 
finite  set,  N  =  8,  the  eigenvalues  for  all  cases  under  consideration 
were  obtained.  A  second  set,  N  =  16,  with  the  Peclet  number,  Pe  =  I , 
was  used  to  check  the  convergence  of  the  method  by  comparing  the  results 
with  the  first  set,  N  =  8. 

The  eigenvalues,  Bn,  were  determined  by  evaluating  the  roots 
of  the  determinant  containing  Bn,  Equation  43,  by  the  Gaussian  reduction 
technique  on  the  maximum  element  of  the  matrix.  Using  an  incremental 
set  of  3  values,  the  determinant  was  evaluated  within  a  selected  Bmax 
range.  The  approximate  location  of  a  root  was  indicated  by  a  change  in 
sign  of  the  value  of  the  determinant  between  two  increments.  The  exact 
values  of  Bn  were  calculated  by  a  controlled  successive  iteration  scheme 
based  on  the  regula-falsi  method.  A  maximum  number  of  twenty  iterations 
were  allowed  with  termination  check  on  eight  digit  accuracy  for  two 
successively  calculated  values  for  each  root  of  the  eigenvalue. 

The  numerical  calculations  were  performed  by  a  CDC  1 40 1 B 
mode!  computer  on  single  precision  level  of  eleven  floating  point  num¬ 
bers.  The  maximum  desired  values  for  the  eigenvalues  were  bounded  by 
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®max  =  2900  for  the  case  of  Pectet  number,  Pe  =  “,  by  Bmax  =  6000  for 
the  case  of  Pe  =  100,  and  by  Bmax  =  600  for  the  cases  of  Pe  =  10  and  I. 
These  ranges  of  Bmax  were  selected  for  calculation  after  the  Initial 
Investigation  of  the  eigenvalues  with  wall  blowing,  RCX  =0.1,  for  the 
case  of  Pe  =  ».  Except  for  the  case,  Pe  =  «,  the  range  of  8max  was 
large  enough  to  yield  8  and  16  roots  for  the  two  finite  sets  of 
trigonometric  functions.  The  calculations  were  carried  out  paramet¬ 
rically  for  Prandt!  number,  Pr  =  0.01,  0.1,  I,  10,  Peclet  number,  Pe  =  «®, 
100,  10,  I,  and  wall  blowing  parameter,  RCX  =0.1,  1,5,  10.  For  a 
given  set  of  parameters  Pr,  Pe,  and  RCX,  the  average  computational  time 
required  for  a  complete  set  of  eigenvalues  was  1.0  minutes  for  the 
function  set  N  =  8  and  4.5  minutes  for  the  function  set  N  =  16. 

The  eigenvalues,  calculated  from  the  determinant.  Equation  43, 
are  presented  in  Tables  I  through  V,  for  parametric  values  of  the  Peclet 
number,  Prandtl  number,  and  the  wall  blowing  parameter,  RCX. 

Inspection  of  the  data  presented  in  Tables  I  through  V  indicates 
that  the  Peclet  number  has  a  strong  effect  upon  the  value  of  the  eigen¬ 
values.  For  the  Peclet  number,  Pe  =  «,  only  three  and  four  eigenvalues 
exist  In  the  range  0  *  Bn  *  2900,  see  Table  I.  Table  II  presents  data 
for  Peclet  number,  Pe  =  100,  and  eight  roots  are  noted  to  exist  In  the 
range  0  -  Bn  -  6000.  For  the  low  Peclet  numbers,  Pe  =  10  and  Pe  =  I , 
the  range  0  -  f$n  -  600  yields  eight  roots  for  the  N  =  8  set  and  sixteen 
roots  for  the  N  =  16  set.  A  comparison  of  a  particular  eigenvalue,  6|, 
for  the  various  Peclet  numbers  shows  that  there  is  a  marked  decrease 
in  the  eigenvalue  magnitude  as  the  Peclet  number  is  reduced.  For  the 
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TABLE  I 

EIGENVALUES  FOR  Pe  =  «.,  N  =  8 


*n 

n 

Pr 

RcX  =  0.1 

>■* 

II 

M  .=  5 

o 

u 

•< 

£ 

1 

0.01 

9.05047(1) 

9.21134(1) 

1.26480(2) 

2.03400(2) 

2 

0.01 

4.316(9(2) 

4.39061(2) 

5.95464(2) 

9.31460(2) 

3 

0.01 

1.02892(3) 

1.04656(3) 

1 .41640(3) 

2.20617(3) 

4 

0.01 

1.88233(3) 

1.91454(3) 

2.58979(3) 

1 

0.  !0 

9.05050(1) 

9.21365(1) 

1.27262(2) 

2.08364(2) 

2 

0.10 

4.31620(2) 

4.39088(2) 

5.96389(2) 

9.37328(2) 

3 

0.10 

1.02892(3) 

1.04659(3) 

1.41739(3) 

2.21245(3) 

4 

0.10 

1.88234(3) 

1.91457(3) 

2.59082(3) 

1 

1.0 

9.05276(1) 

9.44365(1) 

2.04845(2) 

6.81 109(2) 

2 

1.0 

4.31646(2) 

4.41790(2) 

6.88590(2) 

1.51393(3) 

3 

1.0 

(.02895(3) 

1.04947(3) 

» .51625(3) 

2.83737(3) 

4 

1.0 

1.88236(3) 

1.91757(3) 

2.69375(3) 

C 

10.0 

9.27888(1) 

3.22358(2) 

7.07480(3) 

2 

10.0 

4.34301 (2) 

7.09066(2) 

8.43610(3) 

3 

10.0 

1.03170(3) 

1.33630(3) 

1.00050(3) 

4 

10.0 

1.88531(3) 

2.21626(3) 

The  numbers  In  the  parentheses  represent  10^ 
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TABLE  II 

EIGENVALUES  FOR  Pe  =  100,  N  =  8 


n 

Pr 

Rc*  =  o.l 

Rc*  =  1 

Rc*  =  5 

RCA  =  10 

1 

0.01 

8.40486(1) 

8.53422(1) 

1.10910(2) 

1.54939(2) 

2 

0.01 

3.20737(2) 

3.23936(2) 

3.79159(2) 

4.49716(2) 

3 

0.01 

6.05883(2) 

6.10039(2) 

6.78216(2) 

7.58133(2) 

4 

0.01 

9.04658(2) 

9.09314(2) 

9.84295(2) 

1.06941(3) 

5 

0.01 

1.20875(3) 

1.21374(3) 

1.29333(3) 

1.38189(3) 

6 

OcOI 

1.51554(3) 

1.52089(3) 

1.60436(3) 

1.69500(3) 

7 

0.01 

1.83986(3) 

1 .84438(3) 

1.92000(3) 

2.00883(3) 

8 

0.01 

2.15079(3) 

2.15650(3) 

2.24295(3) 

2.33185(3) 

1 

0.10 

8.40488(1) 

8.53621(1) 

1.11522(2) 

1.58011(2) 

2 

0.10 

3.20738(2) 

3.23952(2) 

3.79595(2) 

4.51588(2) 

3 

0.10 

6.05883(2) 

6. 10051(2) 

6.78529(2) 

7.59420(2) 

4 

0.10 

9.04658(2) 

9.09324(2) 

9.84535(2) 

1.07038(3) 

5 

0.10 

1.20875(3) 

i .21375(3) 

1.29352(3) 

1 .38267(3) 

6 

0.10 

i .51554(3) 

1.52089(3) 

1.60452(3) 

1 .69566(3) 

7 

0.10 

1.83986(3) 

1,84439(3) 

1 .92014(3) 

2.00939(3) 

8 

0. 10 

2.15079(3) 

2.15651(3) 

2.24308(3) 

2.33235(3) 

1 

1.0 

8.40684(1) 

8.7345 ! ( 1 ) 

1.69123(2) 

3.92222(2) 

2 

1.0 

3.20753(2) 

3.25542(2) 

4.21879(2) 

6.17795(2) 

3 

1.0 

6.05895(2) 

6.11247(2) 

7.09277(2) 

8.84479(2) 

4 

1.0 

9.04667(2) 

9.10261(2) 

1.00829(3) 

1.16431(3) 

5 

1.0 

1.20876(3) 

1 .21451(3) 

1.31279(3) 

1.45908(3) 

6 

1.0 

1.51555(3) 

1.52154(3) 

1 .62070(3) 

1.75993(3) 

7 

1.0 

1.83987(3) 

1.84995(3) 

1 .93409(3) 

2.06479(3) 

8 

1.0 

2.15079(3) 

2.  55700(3) 

2.25539(3) 

2.38!24(3) 

1 

10.0 

8. 60264 ( L) 

2.60327(2) 

2.12674(3) 

4.72485(3) 

2 

10.0 

3.22333(2) 

4.70391(2) 

2.25070(3) 

4.81060(3) ■ 

3 

10.0 

6.07087(2) 

7.24537(2) 

2.37492(3) 

4.88838(3) 

4 

10.0 

9.05603(2) 

1.00082(3) 

2.50961(3) 

4.96462(3) 

5 

10.0 

1.20952(3) 

1.28920(3) 

2.66530(3) 

5.05057(3) 

6 

10.0 

1 .51619(3) 

1.58482(3) 

2.84497(3) 

5.15268(3) 

7 

10.0 

1.84043(3) 

1.90026(3) 

3.04767(3) 

5.27047(3) 

8 

10.0 

2.15128(3) 

2.20570(3) 

3.27985(3) 

5.42832(3) 

The  numbers  in  the  parentheses  represent  I0^n^ 
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TABLE  III 

EIGENVALUES  FOR  Pe  =  10,  N  =  8 


n 

Pr 

6 

o 

II 

-< 

£ 

RcX  =  1 

RcX  =  5 

V 

II 

o 

1 

0.01 

2.64449(1) 

2.65257(1) 

2.77959(1) 

2.91402(1) 

2 

0.01 

5.84972(1) 

5.85633(1) 

5.95883  0) 

6.06866(1) 

3 

0.01 

9.00428(1) 

9.01076(1) 

9.1  10900) 

9.21679(1) 

4 

0.01 

1 .21497(2) 

1.21562(2) 

1.22564(2) 

1.23616(2) 

5 

0.01 

1.52927(2) 

1.52992(2) 

1.53998(2) 

1.55050(2) 

6 

0.01 

1.84348(2) 

1.84414(2) 

1.85424(2) 

1.86476(2) 

7 

0.01 

2.15787(2) 

2.15852(2) 

2.16852(2) 

2. 17900  (-2) 

8 

0.01 

2.47204(2) 

2.47270(2) 

2.48281(2) 

2.49330(2) 

1 

0.10 

2.64449(1) 

2.65296(0 

2.789350) 

2.95304(1) 

2 

0.10 

5.84972(1) 

5.85652(1) 

5.96375(1) 

6.08832(1) 

3 

0.10 

9.00428(1) 

9.01089(1) 

9.11418(1) 

9.22991(1) 

4 

0.10 

1  .21497(2) 

1.21563(2) 

1.22589(2) 

1.23715(2) 

5 

0.10 

1.52927(2) 

1.52993(2) 

1.54018(2) 

1.55129(2) 

6 

0.10 

1 .84348(2) 

1.84414(2) 

1.85440(2) 

1.86542(2) 

7 

0.  !0 

2,15787(2) 

2.15852(2) 

2. 16866(2) 

2. 17956(2) 

8 

0. 10 

2.47204(2) 

2.47270(2) 

2.48294(2) 

2.49397(2) 

1 

1.0 

2.64488(1) 

2.69155(1) 

3.64735(1) 

5.67200(1) 

2 

1.0 

5.8499)  ( 1 ) 

5.876140) 

6.4373(0) 

7.81454(1) 

3 

1.0 

9.00442(1) 

9.02399(1) 

9.43664(1) 

1.04606(2) 

4 

1.0 

1.21498(2) 

1.21661(2) 

1.25026(2) 

1.33196(2) 

5 

1.0 

1.52928(2) 

1.53072(2) 

1.55975(2) 

1.62815(2) 

6 

1.0 

1.84349(2) 

1.84480(2) 

1.87074(2) 

1.92994(2) 

7 

1.0 

2.15788(2) 

2.15909(2) 

2.18268(2) 

2.23512(2) 

8 

1.0 

2.47205(2) 

2,47320(2) 

2.49522(2) 

2.54255(2) 

1 

10.0 

2.68345(1) 

5.39442(1) 

2.48045(2) 

4.98478(2) 

2 

10.0 

5.86952(1) 

7.59931(1) 

2.54481(2) 

5.01800(2) 

3 

10.0 

9.0(752(1) 

1.02538(2) 

2.64012(2) 

5.06728(2) 

4 

10.0 

1.21596(2) 

1.31139(2) 

2.76692(2) 

5. 13500(2) 

5 

10.0 

1.53006(2) 

1 .60756(2) 

2.92161(2) 

5.22061(2) 

6 

10.0 

1.84414(2) 

1.90931(2) 

3.10020(2) 

5.32329(2) 

7 

10.0 

2.15844(2) 

2.21464(2) 

3.29898(2) 

5.44213(2) 

8 

10.0 

2.47254(2) 

2.52014(2) 

3.51451(2) 

5.57627(2) 

The  numbers  in  the  parentheses  represent  10^ 
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TABLE  IV 

EIGENVALUES  FOR  Pe  =  1,  N  =  8 


en 

n 

Pr 

Rc*  =  0. 1 

■m 

RCX  =  5 

o 

II 

*< 

£ 

1 

0.01 

3.08773(0) 

3.08870(0) 

3.10347(0) 

3.11848(0) 

2 

0.01 

6.23850(0) 

6.23919(0) 

6.24992(0) 

6.26149(0) 

3 

0.01 

9.38181(0) 

9.38248(0) 

9.39287(0) 

9.40383(0) 

4 

0.01 

1.25239(1) 

1.25246(1) 

1.25350(1) 

1.254580) 

5 

0.01 

.1.56658(1) 

1.56665(0 

1.56768(1) 

1.568750) 

6 

0.01 

1.88075(1) 

1.88082  0) 

1.88185(1) 

1.882930) 

7 

0.01 

2. 19492(1) 

2.19499(1) 

2.19026(1) 

2.19709  0) 

8 

0.01 

2.50909(1) 

2.50916(1) 

2.51019(1) 

2.511260) 

I 

0.10 

3.08773(0) 

3.08909(0) 

3.11331(0) 

3. 15762(0) 

2 

0.10 

6.23850(0) 

6.23939(0) 

6.25484(0) 

6.281 15(0) 

3 

0.10 

9.38181(0) 

9.38261(0) 

9.39616(0) 

9.41696(0) 

4 

0.10 

1.25239(1) 

( . 25247 ( E ) 

(.25374(1) 

1,255560) 

5 

0.10 

1.56658(1) 

1.56666(1) 

1.56788(1) 

1.56954(1) 

6 

0.10 

1.88075(1) 

1.88083(1) 

1.88202(1) 

1.883580) 

7 

0.10 

2.E9492([) 

2.19500(1) 

2. 19616(1) 

2. 197650) 

8 

O.iO 

2.50909(1) 

2.509160) 

2.510310) 

2.51 175  0) 

1 

1.0 

3.08813(0) 

3.12823(0) 

3.97665(0) 

5.88148(0) 

2 

1.0 

6,25905(0) 

6.25905(0) 

6.72896(0) 

8.00794(0) 

3 

1.0 

9.38194(0) 

9.39573(0) 

9.71878(0) 

1.064780) 

4 

1.0 

1,25240(1) 

1.25346(1) 

1.278120) 

1, 35038 ( D 

5 

1.0 

1.56659(1) 

(.56744(1) 

1.58745(1) 

1.64640  0) 

6 

1.0 

1,88076(1) 

1.88.48(1) 

1.898360) 

1,9481 !([) 

7 

1.0 

2.19493(1) 

2.19556  0) 

2.210180) 

2.2532)0) 

8 

1.0 

2.50909(1) 

2.50965(1) 

2.52259(1) 

2,56051  0) 

1 

10.0 

3, 32726(0) 

5.85(93(0) 

2.5(580(1) 

5,007480) 

2 

10.0 

6.25836(0) 

7.98581(0) 

2.574410) 

5.037130) 

3 

10.0 

9.39506(0) 

1.06266(1) 

2.668550) 

5,08594(0 

4 

i  0,0 

1.25339(1) 

1.348280) 

2.794910) 

5.15342  0) 

5 

10.0 

1.56738(1) 

1.64431(1) 

2.949400) 

5,238890) 

6 

10.0 

1.88142(1) 

1.94601(1) 

3.12788(1) 

5.34147(1) 

7 

10.0 

2.19549(1) 

2.25511(1) 

3.32649(1) 

5.460220) 

8 

10.0 

2.50959(1) 

2.55841(1) 

3.541850) 

5.59410(1) 

The  numbers  In  the  parentheses  represent  !0'n* 
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TABLE  V 

EIGENVALUES  FOR  Pe  =  1,  N  =  16 


Pn 

n 

Pr 

Rc*  =  o. i 

RcA  s  | 

RcA  s  5 

>* 

ii 

o 

1 

0.0! 

3.08773(0) 

3.08870(0) 

3.10347(0) 

3.11848(0) 

2 

0.01 

6.23850(0) 

6.23919(0) 

6.24992(0) 

6.26149(0) 

3 

0.01 

9.38181(0) 

9.38248(0) 

9.39287(0) 

9.40383(0) 

4 

0.0! 

1.25239(1) 

1.25246(1) 

1.25350(1) 

1.254580) 

5 

0.01 

1.56665(1) 

1.566650) 

1.567680) 

1.568750) 

-6 

0.01 

1.88075(1) 

1.880820) 

1.881850) 

1.882930) 

7 

0.01 

2.19492(1) 

2.19499(1) 

2.19602(1) 

2.19709  0) 

8 

0.01 

2.50909(1) 

2.509150) 

2.51019(1) 

2.511260) 

9 

0.01 

2.82325(1) 

2.823320) 

2.82435(1) 

2.82542  0) 

10 

0.01 

3.13741(1) 

3.13748(1) 

3.13851(1) 

3.139580) 

1 1 

0.0! 

3.45157(1) 

3.451640) 

3,452670) 

3.45374(1) 

12 

0.01 

3.76573(1) 

3.765800) 

3.766830) 

3.767900) 

13 

0.01 

4.07989(1) 

4.079960) 

4.08206(1) 

4.08206  0) 

14 

0.01 

4.39405(1) 

4.394120) 

4.395150) 

4.396220) 

15 

O.Oi 

4.70822(1) 

4.708280) 

4.70931(1) 

4.710380) 

16 

0.01 

5.02238(f) 

5.02244(1) 

5.02347(1) 

5.024540) 

1 

0.10 

3.08773(0) 

3.08909(0) 

3.11331(0) 

3.15762(0) 

2 

0.10 

6.23850(0) 

6.23939(0) 

6.25484(0) 

6.28115(0) 

3 

0.10 

9.38181(0) 

9.38261(0) 

9.39616(0) 

9.31696(0) 

4 

0.10 

1.25239(1) 

1.252470) 

1.253740) 

1.255560) 

5 

0.10 

1.56658(1) 

1.566660) 

1.56788(1) 

1.56954  0) 

6 

0.10 

1.88075(1) 

1.88083(1) 

1.88202  0) 

1.883580) 

7 

0.10 

2.19492(1) 

2.195000) 

2.196160) 

2.19765  0) 

8 

0.10 

2.50909(1) 

2.509160) 

2.510310) 

2.511750) 

9 

0.  JO 

2.82325(1) 

2.823320) 

2.82446(1) 

2.825860) 

10 

O.IO 

3. 13741(1) 

3.13748  0) 

3.138610) 

3.139970) 

II 

0.10 

3.45157  0) 

3.451640) 

3.452760) 

3.454100) 

12 

0.10 

3.76573(1) 

3.765800) 

3.766910) 

3.768230) 

13 

0.10 

4.07989(1) 

4.079960) 

4.08107  0) 

4.082360) 

14 

0.10 

4.39405(f) 

4.394120) 

4.395220) 

4.396500) 

15 

0.10 

4.708220) 

4.708290  ) 

4.709380) 

4.710640) 

16 

0. 10 

5.02238(1) 

5.022450) 

5.023530) 

5.024790) 
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TABLE  V  (Continued) 


Bn 

n 

Pr 

■ 

o 

II 

£ 

Rcx  =  1 

RcA  =  5 

>> 

n 

o 

1 

1.0 

3.08813(0) 

3.12823(0) 

3.97665(0) 

5.88148(0) 

2 

1.0 

6.23870(0) 

6.25905(0) 

6.72896(0) 

8.00794(0) 

3 

1.0 

9.38194(0) 

9.39573(0) 

9.71878(0) 

1.06478(1) 

4 

1.0 

1.25240(1) 

1.25346(1) 

1.27812(1) 

1.35038(1) 

5 

1.0 

1.56659(1) 

1.56744(1) 

1.58745(1) 

1.64640(1) 

6 

1.0 

1.88076(1) 

1.88148(1) 

1.89836(1) 

1.9481 1(1) 

7 

1.0 

2.19493(1) 

2.19556(1) 

2.21018(1) 

2.25321(1) 

8 

1.0 

2.50909(1) 

2.50965(1) 

2.52259(1) 

2.56050(1) 

9 

1.0 

2.82325(1) 

2.82376(1) 

2.83538(1) 

2.86928(1) 

10 

1.0 

3.13742(1) 

3.13788(1) 

3.14844(1) 

3.17912(1) 

II 

1.0 

3.45158(1) 

3.45200(1) 

3.46170(1) 

3.48972(1) 

12 

1.0 

3.76574(1) 

3,76613(1) 

3.77511(1) 

3.80091(1) 

13 

1.0 

4.07990(1) 

4.08027(1) 

4.08864(1) 

4.11255(1) 

14 

1.0 

4.39406(1) 

4.39441(1) 

4.40225(1) 

4.42455(1) 

15 

1.0 

4.70822(1) 

4.70855(1) 

4.71594(1) 

4.73683(1) 

16 

1.0 

5.02238(1 ) 

5.02269(1) 

5.02969(1) 

5.04935(1) 

1 

10.0 

3.12726(0) 

5.85193(0) 

2.51580(1) 

5.00748(1) 

2 

10.0 

6.25836(0) 

7.98581(0) 

2.57441(1) 

5.03713(1) 

3 

10.0 

9.39506(0) 

1.06266(1) 

2.66855(1) 

5.08594(1) 

4 

10.0 

1.25339(1) 

1.34828(1) 

2.79491(1) 

5.15342(1) 

5 

10.0 

i. 56738(1) 

1.64431(1) 

2.94940(1) 

5.23889(1) 

6 

10.0 

1 . 88J42( 1 ) 

1.94601(1) 

3.12788(1) 

5.34147(1) 

7 

10.0 

2. 19549(1) 

2.25111(1) 

3.32648(1) 

5.46021(1) 

8 

10.0 

2.50958(1) 

2.55841(1) 

3.54184(1) 

5.59410(1) 

9 

10.0 

2.82369(1) 

2.86719(1) 

3.77109(1) 

5.74205(1) 

10 

10.0 

3. 1378KI) 

3.17702(1) 

4.01 184(1) 

5.90303(1) 

1  1 

10.0 

3.45194(1) 

3.48762(1) 

4.26215(1) 

6.07600(1) 

12 

10.0 

3.76607(1) 

3.79881(1) 

4.52044(1) 

6.25995(1) 

13 

10.0 

4.08020(1) 

4.11045(1) 

4.78541(1) 

6.45396(1) 

14 

10.0 

4,39434(1) 

4.42245(1) 

5.05602(1) 

6.65715(1) 

15 

10.0 

4.70848(1) 

4.73474(1) 

5.33140(1) 

6.86870(1) 

16 

10.0 

5.02262(1) 

5.04725(1) 

5. 61 084 ( i ) 

7.08786(1) 

The  numbers  in  the  parentheses  represent  I0*n* 
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low  Peclet  number,  Pe  =  I,  the  N  =  16  set  provides  twice  the  number  of 
eigenvalues  as  the  N  =  8  set  In  the  range  0  *  $n  <  600  with  the  approx¬ 
imate  maximum  value,  B^*  =  50.  Since  the  magnitude  of  the  eigenvalues 
Is  the  controlling  factor  In  determining  the  truncation  point  In  the 
Infinite  series  solution.  Equation  26,  the  calculations  for  the  temper¬ 
ature  development  may  be  In  considerable  error.  If  the  number  of  eigen¬ 
values  Is  not  sufficient  to  assure  the  convergence. 

The  wall  blowing  parameter,  RcX,  Is  noted  to  also  have  consider¬ 
able  effect  on  the  magnitude  of  the  eigenvalues  for  a  particular  com¬ 
bination  of  Peclet  and  Prandtl  numbers.  The  magnitude  of  8  tends  to 
increase  as  the  value  of  the  wall  blowing  parameter  Increases.  Inspec¬ 
tion  of  the  data  presented  In  Tables  I  and  II  for  Peclet  numbers,  Pe  =  «* 
and  100,  reveals  that  there  is  a  relatively  large  increase  lr,  the  values 
of  the  eigenvalues  as  the  wall  blowing  parameter  Increases  from  f^X  = 

0.1  to  F^X  =  10.  In  Table  I  for  Prandtl  number,  Pr  *=  I,  the  eigenvalue, 
B|,  Is  seen  to  Increase  from  approximately  90  to  680  as  the  wall  blowing 
parameter  Increases  from  0.1  to  10.  It  is  also  noted  that  the  change  in 
magnitude  of  the  eigenvalue  with  respect  to  the  wall  blowing  parameter  is 
increased  as  the  Prandtl  number  changes  from  Pr  =  0.01  to  Pr  =  10.  These 
same  effects  are  noticed  but  to  a  lesser  extent  In  Tables  Ml  through  V, 
for  Peclet  numbers,  Pe  =  10  and  I.  In  Table  IV  for  Prandtl  number, 

Pr  =  I,  the  eigenvalue,  B|,  Is  seen  to  increase  from  approximately  3  to 
6  as  the  wall  blowing  parameter  Increases  from  0.1  to  10.  For  the 
large  values  of  wall  blowing  parameter  and  Prandtl  number,  I^X  =  10 
and  Pr  *  10,  Tables  III  and  IV  Indicate  a  small  difference  In  value  for 
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successive  eigenvalues  within  the  finite  set  of  functions,  N  =  8.  The 
range  for  Peclet  numbers,  Pe  =  10  and  I,  Is,  respectively,  49.8  <  Bn 
<  55.8  and  50.0  <  8n  <  55.9  for  the  first  eight  eigenvalues. 

One  of  the  interesting  characteristics  of  the  present  method 
is  that  the  eight  eigenvalues,  which  is  the  maximum  number  that  can  be 
evaluated  from  an  eight  by  eight  determinant,  compare  extremely  well 
with  the  first  eight  eigenvalues  evaluated  from  the  N  =  16  set,  for  the 
case  of  Pe  =  I.  Hence,  it  Is  sufficient  to  start  the  solution  for  the 
eigenvalue  problem  with  the  N  =  8  set,  and  later  Increase  the  number  of 
approximating  set  to  N  =  16  If  It  becomes  necessary  to  calculate 
additional  eigenvalues. 
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CHAPTER  V 

ENTRANCE  REGION  TEMPERATURE  DEVELOPMENT 

« 

The  eigenvalues  evaluated  by  the  B.  G.  Galerkln  method  in  the 
previous  chapter  were  used  to  determine  the  temperature  development 
under  constant  wall  temperature  conditions  for  a  specified  uniform 
entrance  temperature,  T0(X2)  =  T*.  The  linear  combination  constants, 
Cj(n>,  associated  with  the  eigenfunction,  Fn(y),  were  evaluated  from  the 
reduced  system  of  linear  equations.  Equation  42.  For  each  particular 
value  of  Bn  a  set  of  Cj^n^'s  was  obtained  by  use  of  the  Gaussian  re¬ 
duction  technique  on  the  pivotal  element  of  the  coefficient  matrix. 

The  analytic  and  approximate  solutions  for  the  fully  developed 
temperature,  ^(y),  are  given  in  Appendix  B.  The  approximate  B.  G. 
Galerkln  method  solution, 

e^ty)  =  PrEc  e!sPrRcXy  ^  Di  sin  (5!) 

1  =  1 


was  used  since  an  extensive  simp J if  I cat ion  in  the  evaluation  of  the 
required  Integrals  for  the  Anfs  was  possible.  In  the  present  study, 
twenty  terms  were  used  to  determine  the  fully  developed  temperature 
given  In  Equation  51.  The  final  form  of  the  dimensionless  excess 
temperature.  Equation  45,  becomes. 


T(y,£)  -  T  T  (y)  -  T 

n<y,s>  =  e<y,e>  -  e„<y>  =  jt—j - - T*  .  T 

'o  'w  'o  ‘w 


(52) 
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T(y,5)  -  Tw 


'o  w 


K  r  N 

y'  An  ,  ei^r^c^y  ^ 

n=l  [  I=| 


Cjiru  sfn  liry  je 


+  PrEc 


Dj  sin  liry 


The  system  of  simultaneous  equations  for  the  evaluation  of  the 


An's  Is  from  Equation  47 


■  ■% 

An  CKm.n)]  = 


Following  the  general  procedure  as  outlined  In  Chapter  III,  the 
necessary  Integrals  for  the  evaluation  of  the  An's  by  the  Welerstrass 
approximation  theorem  are  evaluated  by  analytical  Integration  as, 


1  r  N  /x  IT  n 

I(m,n)  =  f  ^  cj(n'  sin  liry  ^  Cj 

o  L  1  =  1 


sin  Jvy  e^rRc^y  dy 


N  N 


L  Z  CI 

1=1  j=i 


(n)  c(m)  l  !?*£?.  {.,,l+j  ePrRcX  _  , 


r  i  .  X 

I  (PrRcX)2  +  (l-j)V  (Prf^X)2  +  (i+j)V 
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and 


L  M 

L(m)  =J  l-PrEce^MV^  D|  s!n  liry 


N 


^  CJ  sin  Jiry 

■*J-I 


e^rRcXy  dy 


N 

V  c  (ml  - - ~  (-|)J  e&rRcx  -  | 

L.  CJ  (PrRcX)z  t  4jzirz  [  1 

j=l  V 


M  N 


r 


PrEc 


Z  Z; 

1=1  j=l 


(m) 


<-|)<+j  ePrRcX 


I 


iPrl^A)2  +  (I-j)2ir2  (PrRcX)2 


-  I 

■ 


+  <I+J)2ir2  ] 


(55) 


where. 


Cf(m)  =  Cf(n)  =  , 


(56) 


from  the  definition  of  the  eigenfunctions  Rn(y),  and  constants,  An„ 

Since  the  linear  combination  constants,  Cj^n^,  are  completely  determined 
for  the  set  of  eigenfunctions,  the  integrals.  Equation  54  and  Equation 
55,  are  readily  evaluated  for  solution  of  An's  from  the  system  of  linear 
equations. 

The-  system  of  linear  equations  was  solved  by  the  Gaussian  re¬ 
duction  scheme  for  various  values  of  the  parameters  Pr,  PrEc,  and  RCA. 
The  excess  temperature  was  evaluated  at  twenty-one  points  across  the 
channeJ  for  different  axial  locations,  E. 
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The  development  of  the  dimensionless  temperature,  e(y,5),  along 
the  channel  Is  presented  In  Figures  2  through  12  for  various  values  of 
wall  blowing  parameter,  RgX,  Prandtl  number,  Pr,  and  the  product  of  the 
Prandtl  and  the  pseudo  Eckert  numbers,  PrEc.  For  various  axial  loca¬ 
tions,  £,  profiles  are  given  for  values  of  the  Peclet  number,  Pe  =  •, 

100,  10,  and  I. 

The  effect  of  the  wall  blowing  parameter  on  the  dimensionless 
temperature  profiles  Is  shown  In  Figures  2  through  10  for  Prandtl  number, 
Pr  =  I,  The  hydrodynamic  case  with  no  wall  blowing  Is  approximated  by 
the  conditions  presented  In  Figures  2  through  4.  It  is  noted  that  these 
conditions.  In  which  RCA  =0.1,  produce  profiles  which  are  approximately 
symmetrical  about  the  channel  centerline.  Inspection  of  the  data  pre¬ 
sented  In  Figures  5  through  9  shows  that  the  excess  temperature  prof  I les 
become  more  non-symmetrlcal  as  the  wall  blowing  parameter  Increases  from 
=  I  +o  RCA  =  10.  The  maximum  value  of  excess  temperature  for  all 
axial  coordinate  positions,  £,  Is  observed  to  be  displaced  toward  the 
upper  channel  wall  at  y  =  I.  This  result  Is  to  be  expected  due  to  the 
Injection  of  fluid  through  the  lower  channel  wall,  y  =  0,  with  a  temper¬ 
ature  of  Tw.  From  Figure  A-l,  page  73,  the  wall  blowing  parameter  Is 
noted  to  displace  the  velocity  profile  toward  the  upper  channel  wall,  pro¬ 
ducing  greater  velocity  gradients  near  this  wall.  The  displacement  of 
the  velocity  profiles  results  In  greater  heating  due  to  viscous  effects 
In  the  region  near  the  upper  wall  of  the  channel,  which  also  tends  to 
produce  non-*symmetrIca!  temperature  profiles.  Inspection  of  the  data 
Indicates  that  the  wall  blowing  parameter  tends  to  decrease  the  magnitude 
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Fig.  4  Temperature  Profiles;  Pr  =  1,  RCA  <=  0.1,  PrEc  =  1000 
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Fig.  7  Temperature  Profiles;  Pr  =  1,  RCA.  =  1,  PrEc  =  1000 
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Fig.  9  Temperature  Profiles;  Pr  =  1,  ReA  =  10,  PrEc  =  100 
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Fig.  11  Temperature  Profiles;  Pr  =  0.01,  ReA.  =  5,  PrEc  -  100 
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of  the  excess  temperature  profiles  for  a  particular  value  of  5.  This 
can  easily  be  explained  since  Increasing  values  of  RCA  decrease  the  mean 
velocity  of  the  fluid  for  a  given  pressure  drop,  see  Figure  A-2,  page 
74»  which  results  In  a  lesser  effect  due  to  viscous  dissipation. 

Data  are  presented  for  three  values  of  the  product  of  the 
Prandtl  and  pseudo  Eckert  number,  PrEc  =  10,  100,  1000.  This  parameter 
repres_ents  the  ratio  of  the  dissipation  In  the  fluid  to  conduction  due 
to  the  temperature  difference  between  the  entrance  and  the  wall  con¬ 
ditions.  in  each  figure  PrEc  =  10  represents  the  case  of  strong  cool¬ 
ing  effect  due  to  large  temperature  difference  between  the  entrance  and 
the  wall  conditions,  and  PrEc  =  1000  represents  the  case  of  dominant 
overall  dissipation. 

In  the  figures  where  the  wall  blowing  parameter,  R_X  =  0.1  and 
I,  the  excess  temperature  profiles  near  the  entrance  region  exhibit  the 
characteristic  of  a  dip  near  the  channel  centerline.  This  character¬ 
istic  dip  is  the  result  of  the  internal  heat  generation  due  to  viscous 
dissipation  in  the  boundary  layer  region  which  locally  creates  an 
increase  In  the  fluid  temperature.  For  larger  values  of  the  wall  blowing 
parameter,  RCA  =  5  and  10,  this  characteristic  dip  effect  Is  diminished 
due  to  the  displacement  of  the  velocity  profile  toward  the  upper  channel 
wall.  Figure  10,  page  43,  presents  the  case  where  the  major  portion  of 
the  viscous  dissipation  occurs  near  the  upper  wall  of  the  channel. 

Figures  II  and  12,  pages  44  and  45,  are  presented  to  show  the 
effect  of  the  Prandtl  number  on  the  dimensionless  excess  temperature 
profiles.  Comparing  Figure  8,  page  41,  with  Figures  H  and  12  in  which 
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RCX  =  5,  PrEc  =  100,  with  the  Prandtl  number,  Pr  =  0.01,  0.1,  and  I,  the 
temperature  profiles  are  noted  to  become  more  non-symmetrical  as  the 
Prandtl  number  increases  with  the  maximum  temperatures  being  displaced 
toward  the  upper  channel  wall.  The  dimensionless  excess  temperature 
profiles  presented  in  Figures  II  and  12,  pages  44  and  45,  are  noted  to 
be  very  similar;  however,  the  profiles  in  Figure  8,  page  41,  when  com¬ 
pared  to  the  lower  Prandt?  numbers  are  seen  to  be  considerably  differ¬ 
ent.  The  larger  value  of  Prandtl  number,  Pr  =  I,  is  seen  to  have  a 
relatively  strong  influence  on  the  excess  temperature  distribution. 

It  is  interesting  to  note  that  for  the  lower  Pec let  numbers, 

Pe  =  100,  10,  I,  the  complete  temperature  profile  across  the  channel. 

In  a  sense,  lags  compared  to  the  temperature  profiles  at  Peclet  number, 
Pe  =  «.  This  effect  can  easily  be  justified  by  considering  the  fact 
that  the  strong  axial  conduction  at  low  Peclet  numbers  results  in  heat 
transfer  from  the  local  cross  section  to  the  entrance  plane,  in  addition 
to  the  cooling  through  the  channei  walls.  Neglecting  axial  conduction 
by  assuming  the  Peclet  number  infinite  which  is  most  commonly  treated  In 
the  literature  can  introduce  considerable  error  in  the  temperature  pro¬ 
files  for  values  of  the  axial  coordinate,  £  S  1.0. 
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CHAPTER  VI 

CALCULATIONS  FOR  THE  BULK  MEAN  TEMPERATURE 
AND  THE  LOCAL  NUSSELT  NUMBER 


Since  the  temperature  profiles  are  not  necessarily  symmetric 
about  the  channel  centerline,  the  value  of  the  local  Nusselt  number 
based  on  the  heat  transfer  through  the  upper  wall,  y  =  I,  and  through 
the  iower  wall,  y  =  0,  is  expected  to  have  different  values.  The 
local  Nusselt  numbers  for  each  wall  are  defined  as: 


and, 


Nu 


Nu(0,x3)  =  L[Tm(x3)  -  Twr'  <0,x3> 


(L,x3)  =  -L[Tm(x3)  -  VT1  |^<L,x3> 


(57) 


(58) 


where  Tm(x3)  is  the  local  bulk  mean  temperature  of  the  fluid. 


^m^x3^  = 


t  v3‘*2>  T(x2,x3><!x2 
£  Vj(x2)dx2 


(59) 


Equation  59  can  be  non-dimensional i zed  to  yield  a  dimensionless  bulk 
mean  temperature,  6m( C> ,  as 


„  ,,,  _  Tm(X3>  -  Tw  f  u(y)0(y,e)dy 

- -  =  — - 

1  /N  '  1.1  a  I 


o  w 


/'  U(y)dy 
o 


(60) 


The  slopes  of  the  temperature  profile  at  each  wall  become. 
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3T 

9x2  = 

3T  ' 
3x2  (L»x3}  = 

Therefore  the  equations  for  the 


To  "  Tw  30 
L  3y 

TS  ”  Tw  ae 

L  3y 


local  Nusselt 


(0,£) 

(1,0 

numbers  become. 


(61) 

(62) 


(0,?) 

Nu(0,O  =  Nu  =  Nu(0,x3)  = 

|£  (1,0 

Nu(  1 ,0  =  NuU)  =  Nu(I,x3)  =  -  ■ ^(O 


(63) 


(64) 


Hence,  a  modified  form  of  the  local  Nusselt  numbers  In  terms  of  the 
dimensionless  temperatures  becomes, 


Nu(0,O  =^'  u(y)dyj  ^  B  (0,0 


300, 

+  ^(0) 
3y 


/*  u(y)  0(y,Ody  +  f!u(y)  600(y)dy 
o  o 


-I 


(65) 


and 


30w 

+  37-(,) 


Nud.O  =  '  [/'  u<Y>dyj  j^|S.  (1,5) 

f !lo  u(y)n(y,Ody  +  H  u(y)ew(y)dyj 


(66) 


Since  the  general  expressions  for  fi(y,0  and  0OB(y)  are  already 
obtained,  the  terms  in  Equations  65  and  66  can  be  evaluated  for  each  wall 
as: 
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3y-  (0>  =  pf*Ec  Y*  DJ  lir 

1=1 


(72) 


(73) 


An  average  Nusselt  number,  Nu(x3),  based  on  the  net  heat  trans¬ 
fer  from  the  system  can  be  defined  as, 

fI!r<x3)  =  ~AiT 9L-  T  )  (74> 

m  'w' 

where  A  is  the  total  surface  area  of  the  system.  The  average  Nusselt 
number  based  on  the  prevlousiy  defined  Nusselt  numbers  for  each  wail 
becomes , 

mv  =  n^(x3)  =  amm  +2  Nu<L»g>  (75) 

The  dimensionless  bulk  mean  temperature  and  the  local  Nusselt 
number  for  wall  blowing  parameter,  Rc  =0.1,  Prandtl  number,  Pr  =  I, 

Pec  let  number,  Pe  =  *,  and  product  of  the  Prandtl  number  and  pseudo 
Eckert  number,  PrEc  =  0,  are  plotted  in  Figures  13  and  14.  These  con¬ 
ditions  represent  closely  the  hydrodynamic  case  which  neglects  effects  of 
wall  blowing,  viscous  dissipation  and  axial  conduction.  These  data  are 
in  good  agreement  with  those  found  in  Reference  2  which  are  also  plotted 
in  Figures  13  and  14.  The  agreement  of  the  present  data  with  those 
found  in  Reference  2  indicates  confidence  in  the  data  reduction  procedure 
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Fig.  13  Bulk  Temperature  Assuming  No  Viscous  Dissipation  or  Axial  Conduction 
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I  x  id3 

Fig.  14  Local  Nusselt  Number  Assuming  No  Viscous  Dissipation  or  Axial  Conduction 
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used  In  this  study. 

The  development  along  the  axial  direction,  5,  of  the  dimension¬ 
less  bulk  mean  temperature,  0m(?),  and  the  local  Nusselt  numbers, 

Nu(l,£),  and  Nu(0,£),  for  the  upper  and  lower  channel  wall  and  the  aver¬ 
age  Nusselt  number,  fJIT(£),  are  presented  In  Figures  15  through  24  for 
various  values  of  wall  blowing  parameter,  RCA,  Prandtl  number,  Pr, 

Pec  let  number,  Pe,  and  the  product  of  the  Prandtl  number  and  pseudo 
Eckert  number,  PrEc. 

The  effect  of  the  wall  blowing  parameter  on  the  dimensionless 
bulk  mean  temperature  Is  shown  In  Figures  15a  through  18a  for  Prandtl 
number,  Pr  =  I.  It  Is  noted  that  there  Is  little  change  In  the  dimen¬ 
sionless  bulk  mean  temperature  development  as  the  wall  blowing  parameter 
changes  from  0.1  to  I;  however,  the  larger  values  of  wall  blowing  param¬ 
eter,  RCA  =  5  and  10,  tend  to  decrease  the  magnitude  of  the  bulk  mean 
temperature  for  a  particular  value  of  £.  As  in  the  discussion  of  the 
temperature  profiles,  this  effect  can  be  explained  since  Increasing 
values  of  decrease  the  mean  velocity  of  the  fluid  for  a  given  pres¬ 
sure  drop  which  results  in  a  lesser  effect  due  to  viscous  dissipation. 

Inspection  of  the  bulk  mean  temperatures  presented  In  Figures 
17a,  19a,  and  21a  for  wall  blowing  parameter,  RCA  =  5,  and  Prandtl 
number,  Pr  =  0.01,  0.1,  and  I,  Indicates  relatively  small  effect  due  to 
changes  In  the  Prandtl  number.  It  can  be  shown  that  the  fully  developed 
bulk  mean  temperature  for  the  hydrodynamic  case  with  no  wall  blowing  is 

IW  =  0.00476  PrEc  (76) 

'o  "  'w 
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Fig.  15  Heat  Transfer  Parameters;  Pr  =  1,  Rc\  =  0.1 
(a)  Bulk  Mean  Temperature;  (b)  Local  Nusselt  Number 
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(a)  Bulk  Mean  Temperature;  (b)  Local  Nusselt  Number 
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(a)  Bulk  Mean  Temperature;  (b)  Local  Nusselt  Number 


Fig.  19  Heot  Transfer  Parameters;  Pr  =  0.01,  RCA  =  5 
(a)  Bulk  Mean  Temperature;  (b)  Local  Nusselt  Number 
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Fig.  21  Heat  Transfer  Parameters;  Pr  =  0.01,  ReA  =  5 
(a)  Bulk  Mean  Temperature;  (h)  Local  Nusselt  Number 
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(a)  Bulk  Mean  Temperature;  (b)  Local  Nusselt  Number 
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Fig.  23  Heat  Transfer  Parameters;  Pr  =  10,  RCA  -  0.1 
(a)  Bulk  Mean  Temperature;  (b)  Local  Nusselt  Number 
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Fig.  24  Heat  Transfer  Parameters;  Pr  =  10,  RCA  =  1 
(a)  Bulk  Mean  Temperature;  (b)  Local  Nusselt  Number 
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Equation  76  is  in  excellent  agreement  with  the  fully  developed  bulk  mean 
temperatures  found  in  Figure  15a,  page  55,  for  Prandtl  number,  Pr  =  I , 
and  wall  blowing  parameter,  =0.1,  which  approximates  closely  the  no 
wall  blowing  case. 

Figures  15b  through  1 8b,  pages  55  through  58,  with  Prandtl 
number,  Pr  =  I ,  demonstrate  the  effect  of  the  wall  blowing  parameter, 

RqA,  on  the  development  of  the  local  Nusselt  numbers.  As  was  expected 
from  observation  of  the  temperature  profiles,  the  lower  wall  Nusselt 
number,  Nu^*,  Is  decreased  and  the  upper  wall  Nusselt  number,  Nu**^, 
is  increased  as  the  wall  blowing  parameter  increases  from  0.1  to  10. 

From  Figure  18b  in  which  RCA  =  10,  the  upper  wall  Nusselt  number  is 
noted  to  be  approximately  55  while  the  lower  wall  Nusselt  number  Is 
noted  to  be  approximately  1.5  for  the  fully  developed  thermal  conditions. 
Similar  differences  in  magnitude  of  the  Nusselt  number  are  noted  for 
axial  positions,  £,  where  the  thermal  profiles  are  not  fully  developed. 

It  can  be  easily  shown  that  the  hydrodynamic  case  with  no  wall 
blowing  has  an  average  Nusselt  number,  Nu*  =  8.75,  for  fully  developed 
thermal  conditions.  The  data  presented  i.n  Figure  15b  with  Pr  =  I  and 
Rj-A  =  0.1  approximate  the  case  with  no  wall  blowing,  and  the  Nusselt 
number  at  fully  developed  conditions  agrees  reasonably  well  with  this 
value. 

Inspection  of  the  data  presented  in  Figures  15  through  18  indicates 
the  effect  of  the  wall  blowing  parameter  RCA,  on  the  length  of  channel, 

£»  reauired  to  obtain  fully  developed  thermal  conditions.  It  is  seen 
that  increasing  values  of  the  wall  blowing  parameter  tend  to  decrease 
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the  value  of  the  axial  coordinate,  5,  required  to  produce  fully  developed 
bulk  mean  temperatures  and  Nusselt  numbers. 

It  is  noted  from  the  inspection  of  the  curves  for  local  Nusselt 
number  versus  axial  coordinate,  5,  that  a  decrease  In  value  below  that 
for  fully  developed  conditions  exists  for  the  product  of  Prandtl  number 
and  pseudo  Eckert  number,  PrEc  =  10  and  100.  This  dip  effect  in  the 
Nusselt  number  indicates  that  the  temperature  slope  at  the  channel  walls 
is  decreasing  at  a  more  rapid  rate  than  the  decrease  in  the  bulk  mean 
temperature.  This  result  is  due  to  the  large  temperature  difference  between 
the  entrance  plane  and  the  wall  conditions  for  the  lower  values  of  the1 
PrEc  product.  For  small  values  of  axial  coordinate,  5,  there  is  a 
large  temperature  gradient  localized  near  the  channel  walls  due  to  con¬ 
ductive  heat  transfer  to  the  walls.  In  Figure  14,  page  53,  for  the 
case  with  no  viscous  dissipation  effects,  the  Nusselt  number  is  noted 
to  exhibit  a  monoton ical ly  decreasing  characteristic  in  the  thermal 
entrance  region. 

The  effect  of  the  Pec  let  number  on  the  buik  mean  temperature 
development  is  similar  to  that  mentioned  in  the  previous  chapter  con¬ 
cerning  the  dimensionless  temperature  profile  development.  It  is  seen 
that  for  the  low  Peclet  numbers,  Pe  =  10  and  I,  the  attainment  of  fully 
developed  conditions  lags  as  compared  to  the  case  of  Peclet  number,  Pe  = 

"•  Inspection  of  the  data  for  the  local  Nusselt  number  reveals  a  similar 
lag  in  development  for  the  low  Peclet  number  which  again  verifies  the 
importance  of  the  axial  conduction  term  which  Is  usually  neglected  In 
thermal  entrance  region  problems. 
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CHAPTER  VII 
GENERAL  CONCLUSIONS 

From  the  results  of  the  present  study.  It  Is  concluded  that 
mathematically  the  B.  G.  Galerkin  method  has  distinct  computational 
advantages  over  the  classical  analytical  and  numerical  techniques  in 
application  to  the  eigenvalue  problems  associated  with  thermal  entrance 
region  analyses.  The  complete  generality  of  the  method,  pertaining  to 
the  form  of  the  differential  equation  and  the  boundary  conditions,  makes 
It  possible  to  obtain  solutions  for  the  extended  eigenvalue  systems 
resulting  from  the  formulation  of  the  problem  by  considering  the  effects 
of  the  wall  blowing  parameter  and  the  axial  conduction  term.  From  this 
solution,  expressions  for  the  temperature  distribution,  bulk  mean  tem¬ 
perature,  and  local  Nusselt  numbers  have  been  obtained. 

Numerical  results  have  been  presented  to  show  the  effects  of 
the  wall  blowing  parameter,  viscous  dissipation,  and  axial  conduction  on 
the  temperature  and  local  Nusselt  number  development  in  the  entrance 
region  of  a  parallel  plate  channel  under  constant  temperature  wall  con¬ 
ditions.  The  results  verify  the  fact  that  for  small  Peclet  numbers  the 
temperature  profiles  and  local  Nusselt  numbers  deviate  considerably  from 
the  previously  treated  cases  of  infinite  Peclet  number.  The  data  pre¬ 
sented  also  indicate  the  marked  effect  of  the  wall  blowing  parameter  on 
the  temperature  profiles  and  the  heat  transfer  parameters. 
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APPENDIX  A 

SOLUTION  FOR  THE  FULLY  DEVELOPED  VELOCITY  PROFILE 

The  governing  equation  for  the  fully  developed  velocity  pro¬ 
file  is  a  non-homogeneous  linear  differential  equation  as  follows  from 
Equation  13, 

2 

d  u  ,  ,  „  ,  du  ,  . 

^7  (y)  "  Rc*  dy  <y)  =  _l  (Al> 

The  complementary  function  for  Equation  A I  is 

uc(y)  =  c,eRcAy  +  c2  (A2> 

Using  the  method  of  undetermined  coefficients  let  us  find  a  particular 
integral  for  Equation  Al.  We  assume  a  particular  solution  of  the  form 

up(y)  =  c3y  (A3) 

which  yields. 


up(y)  -  <M> 

Application  of  the  boundary  conditions, 

u(0)  =  0,  u(l)  =  0  (A5) 

results  in  the  solution  for  the  fully  developed  velocity  profile  as. 
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"W  =  inr 


' 

y  - 


|  -  eRc*y 

I  -  e^cA 


(A6) 


It  Is  easily  shown  by  application  of  L’HopItal's  rule  that  the  limit  of 
Equation  A6  is. 

Limit  u(y)  =  ^  (y  _  y?)  (A7) 

RcX  -►  0  2 

which  checks  the  ordinary  channel  flow  solution  without  wall  blowing. 
Profiles  of  the  non-dimensional i zed  velocity  as  a  function  of  the  wall 
blowing  parameter,  RqX,  are  presented  in  Figure  Al .  Investigation  of 
the  effect  of  the  wail  blowing  parameter  on  the  fully  developed  velocity 
profile  resulted  In  the  selection  of  f^A  =0.1,  1,5,  and  10  for  further 
calculations. 

The  mean  velocity,  Vm,  is  defined  as 


V”  “  L  t  V3<X2><<X2  *  Vc  u(y)dy 


(A8) 


A  dimensionless  mean  velocity  is  defined  as 


which  becomes 


Ym 

um"^ 


(A9) 


u_  = 
m 


eRcX(2-RcA)  -  RCA  -  2 
2Rc2A2  (I  -  e^) 


<  Af  0) 


A  plot  of  the  non-dimensional ized  mean  velocity  as  a  function  of  the 
wall  blowing  parameter  Is  shown  In  Figure  A2,  page  74. 


72 


AEDC-TR-69-58 


0.13 

0.12 

0.11 

0.10 

0.09 

0.08 
o 
:> 

£j  0.07 

> 

"  0. 06 

0.05 
0.04 
0.03 
0.02 
0.01 
0 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  -0.8  0.9  1.0 

y  -  x2/L 

Fig.  A1  Fully  Developed  Velocity  Profiles  as  a  Function  of  Wall  Blowing  Parameter 
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APPENDIX  B 

SOLUTION  FOR  FULLY  DEVELOPED  DIMENSIONLESS  TEMPERATURE  PROFILE 


The  governing  equation  for  the  fully  developed  temperature 
profile  is  a  non-homogeneous  linear  differential  equation  as  follows 
from  Equation  20, 

d2e  de  r  h  i .  l2 

<y>-PrRcA  (y)  =  -PrEc[^  (y)J  (Bl ) 

with  boundary  condition^ 


QJ\)  =  0  ,  0W(O)  =  0 


(B2) 


Since  the  analytic  solution  for  the  fully  developed  tempera¬ 
ture  profile  is  a  very  complicated  expression,  an  extensive  simplifica¬ 
tion  in  evaluation  of  the  integral,  Equation  49,  is  possible  by  using 
the  G.  B.  Galerkin  method  to  obtain  a  series  solution  for  Equation  Bl. 
Defining  a  new  function  as 


GJy> 


which  is  substituted  Into 


d2G 


Pr^Rr 


dy 


2=<y>-—T 


=  e-J5prRcXy  Bjy) 

Equation  Bl  to  yield, 

2*2 

- Gw(y)  =  -PREce"3SPrRc^y 


(B3) 


(B4) 
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Assume  an  approximate  solution  of  Equation  B4  as 


M 


G.(y)  =  ^  Dj"  sin  liry 
1  =  1 

Substitute  Equation  B5  Into  Equation  B4  which  gives. 


M 

z 

1=1 


r 


Df  K 


_  pr2R  2»2 

-l2ir2  sin  liry  -  - _£ - sin  liry 


)■ 


(B5) 


PrEc  e'*yarRcA 


du  ,  . 
Hy  (Y} 


1 2 


(B6) 


Following  the  method  of  B.  G.  Galerkin  hence  multiplying  Equation  B6  by 
♦^Cy)  =  sin  kiry  and  integrating  over  the  range  0  ^  y  ^  I, 


M 


z  »r 


i  =  l 


-I2*2  - 


Pr2^2*2 


I,(k,i)  =  -PrEcI2(k) 


k=  1,2,3,  .o,,  M 


(B7) 


where  the  integrals  are  given  by: 


Ii(k,|)=  J  siniirysin  kiry  dy  =  —  6.. 

1  o  7  2  ik 


(B8) 


76 


AEDC.TR. 69.58 


I,(ki  =  /' 


(y) 

dy 


sin  kiry  dy  = 


4it 

W 


2RrX 


Pr2Rc2X2+4k2ir2 


|_l-eRc*J  L(2RcX-PrRcX)2+4k2Tr2 


I  _  e"^rRcX  (-I)1 


|_e^RC^_y>rRCX^  (-1)^ 


Rr2X2 


a-/0*)2 


_(4RcX-PrRcX)2+4k2ir2 


l_(_  I )  ke(2RcX-JsPrRcX) 


(B9) 


Equation  B7  Is  a  non-homogeneous  system  of  linear  equations  which  can 
be  solved  to  determine  the  Dps.  Since  I,(k,l>  is  the  Kronecker  delta 
Equation  B7  can  be  written  as  follows. 


Di 


-2PrEc  I2U) 


2  2  pr2Rc2X2 


- 


if  Df  is  redefined  as. 


°'T 


-2I2(i) 


Pr2 R  2*2 


then. 


(BIO) 


(Bl  I  ) 


GJy)  =  PrEc 


I 

t  _  • 


Dj  sin  i Try 


(BI2) 
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Fron  Eouation  H3  the  fullv  developed  dimensionless  temperature  profile 
becomes 

».<y>  -  eVri,cXy  G.(yj  ■ 


The  approximate  solution  for  the  present  study  was  obtained  by  using 
twenty  terms.  M  =  20,  in  the  series  solution. 
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